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We calculate the Z2 topological invariant characterizing a one dimensional topological supercon- 
ductor using a Wess-Zumino-Witten dimensional extension. The invariant is formulated in terms 
of the single particle Green's function which allows to classify interacting systems. Employing a 
recently proposed generalized Berry curvature method, the topological invariant is represented in- 
dependent of the extra dimension requiring only the single particle Green's function at zero frequency 
of the interacting system. Furthermore, a modified twisted boundary conditions approach is used 
to rigorously formulate the topological invariant for disordered systems. 



Introduction. The key feature of a topological 
superconductor (TSC) in one spatial dimension (ID) 
[T}]3] is a holographic single Majorana bound state 
(MBS) associated with each of its ends which is topolog- 
ically protected by particle hole symmtery (PHS). The 
recently proposed realization of the ID TSC in an InSb 
nanowire with strong spin orbit interaction (SOI) and 
proximity induced s-wave superconductivity [H |3] is so 
far the most promising candidate as to its experimental 
feasibility. Very recently, first experimental signatures of 
MBS have been reported by several experimental groups 
[5H7], however, it should be mentioned that alternative 
explanations for robust zero bias resonances not owing 
to Majorana zeromodes have been brought forward 
[Si [9] . Note that PHS in this class of systems does not 
imply a true limitation on the band structure of the 
normal- conducting spin orbit coupled quantum wire but 
is emergent from the Bogoliubov-deGennes (BdG) mean 
field description of superconductivity. An additional 
chiral symmetry present in the ideal model systems 
proposed in Refs. [2j promotes the Z2 invariant 
characterizing the presence of an unpaired MBS to a 
Z invariant [101 111) counting the number of zeromodes 
at each edge. However, perturbations modifying the 
SOI as well as magnetic impurities can break the chiral 
symmetry and gap out paired MBS. The influence of 
interactions on the topological classification of chiral 
ID systems has been analyzed from a matrix product 
state perspective [T^]. Very recently, the classification 
of a chiral ID system in terms of its single particle 
Green's function has been reported jlS, . The robustness 
of the TSC phase in interacting nanowircs has been 
investigated using renormalization group (RG) methods 

[HHS]. 

In this work, we go beyond the Z2 classification in 
terms of the Pfaffian of the Hamiltonian in Majorana 
representation reported in Ref. [T] by including disorder 
and adiabatic interactions into the classification scheme. 
We present the Z2-classification of the ID TSC phase 
which, without additional symmetries, generically 
belongs to the Cartan-Altland-Zirnbauer class D [T7] . 



in terms of its single particle Green's function. In a 
first step, we work out explicitly a dimensional exten- 
sion procedure for a realization of the ID TSC in a 
nanowire with strong SOI and proximity induced s-wave 
superconductivity [21 [3]. This allows us to reduce the 
topological classification of the ID TSC to the quantum 
anomalous Hall (QAH) effect [T5] of the extended 2D 
system. This procedure fits into the general classification 
framework of topological field theory (TFT) proposed 
for time reversal invariant topological insulators in Refs. 
(T9| [20] and allows for a reformulation of the invariant 
in terms of the single particle Green's function. Upon 
switching on interactions adiabatically, our classification 
remains valid for Luttinger liquid like interactions as 
argued in Ref. [21]. Thereafter, we employ a recently 
proposed generalized Berry curvature method |22j 
showing that the interacting invariant can be expressed 
in terms of the zero frequency single particle Green's 
function of the physical ID system, which is indepen- 
dent of the previously introduced extra dimension. 
Finally, we demonstrate how a hybrid approach of 
twisted boundary conditions (TBC) in the physical 
dimension and periodic boundary conditions in the 
extra dimension can be used to additionally include 
disorder at the level of the bulk topological invariant, 
i.e. without probing the presence of unpaired MBS, 
quantized zero bias resonances, or other finite size effects. 

Model of the ID TSC. A lattice model of the ID TSC 
[H|TT] can be cast into the form H = J 'i/^HBdG'^, where 
the basis is chosen such that ^I' = {ip^, V^, V'j, V'|)- In this 
basis 

HBdG = (Jt Jjj?j ■ 

For the ID TSC, Hoik) = Cfe + Ba^ + usm{k)ay is the 
Hamiltonian of a single channel quantum wire in the pres- 
ence of a i?-field induced Zeeman splitting and Rashba 
SOI. The proximity induced s-wave superconducting gap 

is of the form S — and ^/t = 1 — cos k — fi. 

Introducing the set of Pauli-niatrices for the particle 



hole pseudo spin, the BdG Bloch Hamiltonian reads 

HBdcik) = {(.k + Ba.j; + usm{k)ay) + AayTy. (1) 

In this representation, the PHS operation has the intu- 
itive form 

C = T,K, (2) 

where K denotes complex conjugation. Let us very 
briefly review the salient physics starting from the 
continuum model obtained from Eq. ([l]) by substituting 
sm{k) k, cos{k) i-^ 1 - ^. For B = A = 7^ w, the 
band structure consists of two particle hole symmetric 
copies (emergent from the BdG picture) of the shifted 
Rashba parabulae. The lattice regularization in Eq. 
^ is introduced to make the topological invariants 
well defined. A ^ gaps out the system in its entire 
Brillouin zone (BZ). For small k this gap competes with 
a Zeeman gap due to _B 7^ leading to a band inversion 
at = ^i^ + A^. For > + A^ we have a TSC with 
a single MBS associated with each end of a finite wire. 

Dimensional extension. Following the general out- 
line in Ref. |19j . we explicitly perform a dimensional ex- 
tension introducing an extra coordinate v, thus reducing 
the topological classification of the noninteracting model 
to the analysis of the QAH effect of the extended 2D 
system. The idea is quite simple: Our system cannot 
be deformed into a trivial ID insulator without break- 
ing PHS which provides the topological protection of the 
TSC phase. However, breaking this symmetry we can 
deform the TSC, say upon varying v from to tt, into a 
trivial ID superconductor without ever closing the bulk 
gap of the instantaneous system. If we were to investi- 
gate the QAH features, i.e. the first Chern number C, of 
the 2D system resulting from such an interpolation, the 
integer would totally depend on our choice of the inter- 
polation. However, performing the interpolation also for 
V e [— TT, 0] in such a way that the resulting 2D system 
is 27r-periodic in v, PHS preserving, and insulating, C is 
well defined up to even integers. 

This means that a Z2 information 1^ = C(mod 2) is well 
defined and only depends on the physical ID system. It is 
worth noting that finding a suitable interpolation is non- 
trivial and requires some insight into the physical mecha- 
nisms underlying the model. In the following, we will ex- 
plicitly present an extension which works for a generic ID 
TSC and can also be used later on for the disordered in- 
teracting system. The guiding physical idea is as follows: 
Switch on a particle hole breaking gap ^ sin(w)T2, which 
will keep the gap open for v 0,tt (mod 2tt). Destroy 
the band inversion by a term ~ /3(2 — 2 cos v)ayTy which 
vanishes for the physical model (v = 0) and will for suf- 
ficiently large /3 produce a trivial superconducting phase 
for V = ±7r where it enhances the superconducting gap 
by 4/3. In summary, the Wess-Zumino-Witten (WZW) 
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FIG. 1. (Color online) 3D plot of the Berry curvature 
J-{k,v) for nontrivial parameters {B — 1.5) left and trivial 
parameters {B — 0.4) right, u = A — 2/3 — 1, ^ = in both 
plots. 

pi] extended Hamiltonian reads 

U-wzwik^v) = i?BdG(fc) + sin(u)rj. + /3(2 - 2 cos v)ayTy. 

(3) 

Integrating the Berry curvature T of this Hamiltonian 
over the {k,v) BZ indeed yields 

C^^f F = e{B^-A^-i?), (4) 
Jbz 

valid for parameters close enough to the band inversion 
that no artificial level crossings which depend on the 
details of the lattice regularization occur. 

In Fig. [1] we compare the Berry curvature of an 
extension of a nontrivial ID TSC with that of a trivial 
superconducting wire. In the extra dimension w, the 
modulus of the curvature is smoothly decaying without 
any notable difference between the trivial and the 
nontrivial case. This is reflected in our derivation 
below, which shows that the topological invariant of the 
translation-invariant system can be defined in terms of 
its single particle Green's function without reference to 
the extra dimension. Note that this picture changes in 
the framework of TBC as introduced to account for the 
presence of disorder as is discussed below and shown in 
Fig. [I 

Single particle Green's function topology. We now 
discuss the possiblity to include interactions into the 
proposed classification scheme by generalizing the Chern 
number of the noninteracting system to a topological in- 
variant of the single particle Green's function in com- 
bined frequency momentum space, as has been proposed 
for time reversal invariant topological insulators in 2D 
and 3D |20j . In even spatial dimension 2n, it has been 
shown p5H27] that the homotopy of the single parti- 
cle Green's function determines the Hall conductivity 
(^xy (ji = 1) and its higher dimensional analogues, re- 
spectively, of an insulating fermionic system. It has been 
explicitly demonstrated that this representation of 
Gxy adiabatically connects to the nth-Chern number as- 
sociated with the Berry curvature of a Bloch Hamiltonian 
in the non-interacting limit. The adiabatic assumption 
in this case implies that the gapped ground state of the 
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interacting system can be continuously connected to the 
manifold of occupied single particle levels of the nonin- 
teracting system which is also separated from the empty 
states by a finite energy gap. The most prominent exam- 
ple of interactions which are non-adiabatic in the sense 
just defined lead to the ^ fractional quantum Hall (FQH) 
effect in a partially filled lowest Landau level (LLL). In 
this case, the noninteracting system is gapless and the 
single particle states that span the LLL split into v de- 
generate ground states when TBC are applied [33]. This 
ground state degeneracy is at the basis of the concept 
of topological order j^Sj which rigorously classifies FQH 
systems. 

In ID, interactions play a peculiar role generically lead- 
ing to non-Fermi-liquid behavior. From a viewpoint of 
perturbation theory interactions are therefore considered 
to be non-adiabatic in ID, as no meaningful quasiparti- 
cles can be defined. However, it has been argued [21] 
that the Fermi surface properties as described by the 
momentum space topology of the single particle Green's 
function are still adiabatically connected to those of the 
non-interacting system. 

Rewritten in terms of the single particle Green's func- 
tion G{uj,p) of the extended system, i.e. Gq{uj,p) = 
{iuj — H-wzwip)) ^ for the special case of the non- 
interacting system, the Z2 invariant v reads 



247r2 



where G, 



[ Tr[GG-^GG-^GG-^] (mod 2), (5) 



d^G-\ 



- "M-^ ' 1^ = = i^,k,v. Very 

recently, it has been generally shown [55] that for 
an interacting system, an invariant of this form can 
be simplified by introducing a generalized Berry cur- 
vature T = —iJ^B.-zerosi'^iP^'^l) ^ d,\p,a) associatcd 
with the fictitious non-interacting Hamiltonian H = 
—G~^{0,p), which takes into account the eigenvectors 
\p, a) of G~^{Q,p) with positive eigenvalues, the so called 
R-zeros [22]. The Z2 invariant then takes the form of a 
generalized Chern number, i.e. 



1 

2n 



P (mod 2). 



(6) 



BZ 



As has been shown for the noninteracting case in Ref. 
[TO] the Z2 classification of the particle hole symmetric 
ID system can then be further simplified to 



2P(0)(mod 2) 



27r 



dkA{k) (mod 2), (7) 



where P(0) is the charge polarization of the physical 

ID system and A{k) — — * X]fl-zeros(^' '^Ml'^' o^) 
generalized Berry connection restricted to the physical 
system at u = 0, i.e. at p = (fc,0). Note that this 
general form does no longer depend on the dimensional 
extension procedure and can be calculated once the 



zero frequency single particle Green's function 0(0, k) is 
known. Finally, the Z2-invariant can be practically 
calculated by formal analogy to the non-interacting 
case by calculating the Majorana number [1] defined in 
terms of the Pfaffian of the fictitious non-interacting 
Hamiltonian H in Majorana representation. 

Disorder and Twisted Boundary Conditions. Our for- 
mulation so far has been relying on translation-invariance 
which implies the existence of a BZ. This description will 
thus no longer be applicable in the presence of disor- 
der. To this end, the concept of TBC has been intro- 
duced to topologically classify quantum Hall systems in 
the absence of translation-invariance [23]. As long as 
a bulk mobility gap is present, the Green's function is 
exponentially bounded in real space for energies in this 
gap. Under these conditions, Niu et al. [23] show that 
the Hall conductivity can be represented as a constant 
ground state Berry curvature with the wave vector re- 
placed by the twisting angles 6*, 4> of the TBC. In this 
formalism, the Hall conductivity axy reads 



'xy 



27riGo 



■ de^ 86' 



1 8^)0 dtpo , 
' 86^ 89 ' 



2'KiGoJ^e<h ~ Gc 



iFa 



rp2 



2Tr ' 



(8) 



where ipo denotes the ground state wave function, Gq = 

2 1 

^ = 2^ is the quantum of conductance. In last equal- 
ity of Eq. ([s]), the independence of J^O(t> on the twisting 
angles [23] has been used to make the topological quanti- 
zation of axy manifest by representing it as Go times the 
Chern number of the f7(l)-bundle over the torus of 
the twisting angles {9, 4>). Since in this work, we consider 
a disordered ID system, we can without loss of general- 
ity assume translation-invariance in the extra dimension. 
Integrating over the momentum v associated with the 
direction of translational invariance is equivalent to eval- 
uating Eq. Q for = for the special case of a system 
with translational invariance in a;-direction. Within this 
hybrid approach of twisted boundary conditions in the 
physical dimension and w-momentum integration in the 
extra dimension, the Chern number of the extended 2D 
system in the presence of this stripe-like disorder can be 
expressed as 



C = 



dv 



d<P 



dv F^A 



(9) 



where T is the Berry curvature on the "mixed" torus 
defined by the wave vector in v-direction and the twist- 
ing angle (f) of the TBC imposed in the physical direction. 
The first equality sign in Eq. (|9| makes the integer quan- 
tization of our topological invariant manifest, whereas 
the second equality sign follows from the independence 
of J^^ dvPyrfi = -^Btplg^Q of the twisting angle (j). The 
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FIG. 2. (Color online) Plot of the Berry curvature 
J, , for nontrivial parameters (A = 0.7, -y — u — B = 

i<#>=0 

P = 1, fi — 0) (purple) and trivial parameters ( A = 0.7, 7 = 
5, u = B = /3 = 1, /I = 0) (blue). Wire length 100 sites in 
both plots. 



main advantage of Eq. (|9| as compared to the general 
2D case (see Eq. ([8|) is that only the cigenstates of a ID 
system have to be calculated to evaluate the topological 
invariant which is numerically less costly. This program 
allows for a topological classification of disordered sys- 
tems with periodic boundaries, i.e. without explicitly 
probing the presence of unpaired MBS. The influence of 
disorder on the presence of MBS has been studied using 
a scattering matrix approach [29l |30] . 

In Fig. [2] we show the mixed Berry curvature 



b=0 



for a weakly disordered system (7 = 1) in the 



topologically nontrivial phase and a strongly disordered 
trivial system (7 = 5). Here, 7 is the strength of a scalar 
Gaussian onsite potential. Note that in contrast to the 
translation-invariant case (see Eq. Q) the topology is 
determined by the w-dependence of the mixed Berry cur- 
vature. Finally, we would like to point out that even for 
7 = 1, the onsite potential fluctuations significantly ex- 
ceed the bulk insulating gap of the ID TSC. The disorder- 
induced transition from nontrivial to trivial takes place 
at disorder strengths which are, depending on the other 
model parameters typically three to five times larger than 
the bulk gap which is in agreement with recent results 
obtained from level spectroscopy in ID TSC with closed 
boundary conditions |31j . 

Conclusions. We have constructed a dimensional 
extension which reduces the topological classification 
of the ID TSC phase to the QAH effect in two dimen- 
sions. This approach is ready-made to rephrase the 
invariant in terms of the single particle Green's function 
of the extended system. Using a generalized Berry 
curvature method, the invariant can be simplified to 
the noninteracting classification scheme with a fictitious 
non-interacting Hamiltonian defined in terms of the 
Green's function of the physical ID system at zero 
frequency. To obtain a well defined topological invariant 
for disordered systems, the method of TBC has been 
modified to a hybrid approach accounting for the partial 



translation invariance in the extra dimension. 
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